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ABSTRACT 

We report the discovery of a unique gravitational lens system, SDSS J2222+2745, producing five 
spectroscopically confirmed images of a z s — 2.82 quasar lensed by a foreground galaxy cluster at 
Zi = 0.49. We also present photometric and spectroscopic evidence for a sixth lensed image of the 
same quasar. The maximum separation between the quasar images is 15" 1. Both the large image 
separations and the high image multiplicity of the lensed quasar are in themselves rare, and observing 
the combination of these two factors is an exceptionally unlikely occurrence in present datasets. This 
is only the third known case of a quasar lensed by a cluster, and the only one with six images. 
The lens system was discovered in the course of the Sloan Giant Arcs Survey, in which we identify 
candidate lenses in the Sloan Digital Sky Survey and target these for follow up and verification with 
the 2.56m Nordic Optical Telescope. Multi-band photometry obtained over multiple epochs from 
September 2011 to September 2012 reveals significant variability at the ~ 10 — 30% level in some of 
the quasar images, indicating that measurements of the relative time delay between quasar images 
will be feasible. In this lens system we also identify a bright (g = 21.5) giant arc corresponding 
to a strongly lensed background galaxy at z s — 2.30. We fit parametric models of the lens system, 
constrained by the redshift and positions of the quasar images and the redshift and position of the 
giant arc. The predicted time delays between different pairs of quasar images range from ~ 100 days 
to ~ 6 years. 

Subject headings: galaxies: clusters: general - - gravitational lensing - - quasars: individual 
(SDSSJ2222+2745) 



1. INTRODUCTION 

Strongly gravitationally lensed quasars with image sep- 
arations > 10" are exceptionally rare. Such large im- 
age splittings require a lensing mass corresponding to a 
galaxy cluster, and a fortuitous close line-of-sight align- 
ment of such a mass distribution and a background 
quasar. In fact, only two such systems have previ- 
ously been discovered, SDSS J1004+4112 (Inada et al. 
2003) with a maximum image separation of 14'/ 6 and 
SDSS J1029+2623 (Inada et al. 2006) with a maximum 
image separation of 22'.'5. Such systems provide unique 
information on the abundances and internal structure of 
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galaxy clusters, yielding stringent tests of the current 
ACDM paradigm, both through their observed abun- 
dances (e.g., Hennawi et al. 2007), the relative fractions 
of different image multiplicities (e.g., Oguri & Keeton 
2004), and from detailed models of the cluster mass dis- 
tribution based on the available constraints (e.g., Sharon 
et al. 2005; Oguri 2010; Oguri et al. 2012), including mea- 
sured time delays between the quasar images (Fohlmcis- 
ter et al. 2008, 2012). Being effectively point sources, the 
positions of multiply lensed quasar images provide mass 
modelling constraints which are free from the ambiguities 
which sometimes plague the identification of multiple im- 
ages of strongly lensed galaxies, particularly the fainter 
de-magnified images located closest to the mass center. 

As shown by Refsdal (1964), the relative time delays of 
multiple quasar images visible around the lensing mass 
center can provide an estimate of the Hubble parameter 
Hq. However, as Hq is presently increasingly well con- 
strained by a host of methods, including gravitational 
lensing time delays from galaxy-scale lenses, the main 
interest of quasar time delays are related to their ability 
to provide unique detailed constraints on the slope of the 
projected matter density distribution over the range of 
radii where the quasar images are found (e.g., Refsdal 
2004). 

In this paper, we present the discovery of a third lens- 
ing system where a quasar is being multiply lensed by a 
galaxy cluster, but with a higher image multiplicity than 
the previously known such systems. The discussion in 
subsequent sections of this paper describes the evidence 



2 



Dahle et al. 



for the lensing interpretation of this system and the iden- 
tification of an increasing number of components in the 
lensing system, through analysis of progressively deeper 
photometric and spectroscopic data, coupled with para- 
metric modelling of the lens. 

In § [2] we describe the initial identification of 
SDSS J2222+2745 as a likely gravitational lens system 
in the course of a larger survey for strong cluster lenses. 
The spectroscopic confirmation of the t hree brigh test 
lensed images of the quasar is outlined in § I3.1[ and § |3.2| 
details the evidence for three additional fainter quasar 
images from deeper imaging of the field, coupled with 
further spectroscopic observations. In § |3.3| we report 
photometry of the quasar images, some showing signifi- 
cant brightness changes over a time span of one year. In 
§ |3.4| we provide a summary of the observed properties 

01 the lensing cluster, and in §[4] we develop a parametric 
model of the lens system, which is used to predict time 
delays. Our main results are summarized in § [5] Un- 
less otherwise noted, throughout this work we assume an 
TO = 0.3, £1a = 0.7, H = 70 km s -1 cosmology. 

2. LENS CANDIDATE SELECTION AND INITIAL IMAGING 

We have undertaken a search for cluster lenses, the 
Sloan Giant Arcs Survey (SGAS; Hennawi et al. 2008; 
Bayliss et al. 2011; Gladders et al. in prep. 2012), based 
on public data from the Sloan Digital Sky Survey (SDSS; 
York et al. 2000). Potential lines of sight for lensing 
by clusters are identified by running an automated clus- 
ter finding algorithm on the SDSS imaging data. This 
is followed by visual inspection, by multiple persons in 
our collaboration, of color images of these lines of sight, 
generated from the SDSS data. New lens candidates 
are flagged and given a grade according to the likeli- 
hood of being a true lens, and a final grade is calcu- 
lated for each lens candidate. The list of lens candi- 
dates is then checked against public data bases such as 
the NASA/IPAC Extragalactic Database, and remaining 
candidates after removal of known lenses are targeted for 
deeper follow-up imaging. In September 2011, we tar- 
geted such lens candidates selected from the SDSS Data 
Release 8 (Aihara et al. 2011) using the MOsaic CAm- 
era (MOSCA) at the 2.56m Nordic Optical Telescope 
(NOT). MOSCA is a 2 x 2 mosaic of four 2048 x 2048 
Loral-Lesser thinned CCDs with excellent blue sensi- 
tivity (Quantum Efficiency > 90% across the SDSS g- 
band), mounted at the Cassegrain focus of the telescope. 
MOSCA is normally used in 2 x 2 binned mode, yielding 
a pixel scale of 0.217"pixel _1 and a total field of view of 
7:7 x 7:7. 

SDSS J2222+2745 was imaged using MOSCA in 
September 2011 (see Table fl] for details) with our nomi- 
nal survey exposure time of 2 x 300s in the g-band. Based 
on a visual inspection of the raw data, it was noted as a 
probable lens, and subject to immediate follow-up with 

2 x 150s exposures in each of the SDSS r- and i-bands. 
In addition to the apparent blue arc south of the clus- 
ter core, three blue stellar images (labelled "A","B", and 
" C" in Figure [I]) located at similar cluster-centric radii 
as the arc, were considered as candidate multiple images 
of a gravitationally lensed quasar. In addition, three red 
galaxies, labelled "G1-G3" in Figure [TJ were seen in the 
cluster core. Publicly available photometric redshift es- 
timates from SDSS (see Table [5| indicated that all three 



galaxies are consistent with being located in a cluster at 
2-0.50. 

The stellar source (labelled "WD" in Figure []} of mag- 
nitude g = 21.09, similar in apparent brightness to A and 
B, seen superposed on the giant arc, was also noted as 
a candidate quasar image. Although this image has a 
significantly redder color (g — i = 1.33) than A and B, it 
might be considerably reddened by extinction in galaxies 
(the galaxy corresponding to the arc being one potential 
such candidate) along its line of sight, and it was kept as 
a target for future spectroscopic follow-up. 

3. FOLLOW-UP OBSERVATIONS AND RESULTS 

3.1. Spectroscopic confirmation of the three brightest 
quasar images 

The spectroscopic observations described in this paper 
were obtained using the Andalucia Faint Object Spec- 
trograph and Camera (ALFOSC) at the NOT. We used 
grism # 4, covering a wavelength range 3800A - 9000A 
with a resolution of about R = 350 and dispersion 3.0 A 
pixel -1 . The ALFOSC pixels were binned by a factor of 
two in the dispersion direction but kept unbinned at the 
native pixel scale of 0'.' 188 pixel -1 in the spatial direction. 
The slit width was l'/3. Redwards of - 7000 A strong 
fringing is present in the spectra, effectively restricting 
the upper limit of the spectral range to this value unless 
special care is taken during the observations to suppress 
the effects of fringing. An order-sorter filter was used to 
suppress second-order contamination in the spectra. 

On August 21-22, 2012, a series of ALFOSC spectra 
were obtained in the area of SDSS J2222+2745: Three 
1200s spectra were obtained with the rotation angle set 
to provide simultaneous coverage of the A and B im- 
ages. In addition, two single spectra of 1000s exposure 
time were obtained to simultaneously cover image C and 
galaxy Gl and image B and image C, respectively. Fi- 
nally, the blue giant arc and the stellar source super- 
posed on it (labelled "WD" in Figure fl]) were observed 
for a total exposure time 3 x 1000s. The spectra were 
wavelength calibrated using He+Ne arc lamps and flux 
calibrated based on observations of the spectrophotomet- 
ry standard star BD+33 2642. As shown in Figure [2] the 
three images A,B, and C all have identical spectra mod- 
ulo noise, containing typical spectral features of a quasar 
(e.g., Vanden Berk et al. 2001) with a set of prominent 
emission lines, redshifted to z = 2.82. This demonstrates 
that the lensing interpretation is correct for these three 
quasar images. 

Despite having a significantly redder color, the stel- 
lar source superposed on the giant arc was also targeted 
for ALFOSC spectroscopy as a potential additional im- 
age of the quasar. Although a significant continuum was 
detected, emission lines were completely absent at the 
wavelengths of the prominent quasar emission lines in the 
A,B, and C component spectra. A more careful analysis 
of the spectrum reveals the typical signatures of a DZ 
white dwarf (e.g., Koester et al. 2011). 

3.2. Identification of the D, E and F images 

Having confirmed the three images A-C as gravitation- 
ally lensed images of the same quasar, two alternative im- 
age configurations may be considered: SDSS J2222+2745 
is either a naked cusp lens with three quasar images, 
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Table 1 

Log of NOT/MOSCA imaging observations 



Epoch 


Filter 


Exposure time 


FWHM 




2011 Sep. 24.93 UT 


9 


2 X 300s 


0''87 


2011 Sep. 24.94 UT 


r 


2 X 150s 


0'.'81 


2011 Sep. 24.95 UT 


i 


2 x 150s 


0'.'74 


2012 Sep. 13.97 UT 


u 


3 x 600s 


CK'81 


2012 Sep. 12.96 UT 
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3 x 600s 


0'.'70 


2012 Sep. 13.04 UT 
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3 x 300s 


0'.'68 


2012 Sep. 12.95 UT 
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3 x 300s 


0'.'61 


2012 Sep. 15.98 UT 
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Figure 1. NOT/MOSCA gri color composite image of SDSS J2222+2745 from our initial epoch of photometry (see Table 



similar to SDSS J1029+2623 (Oguri et al. 2008) or a 
five-image system similar to SDSS J1004+4112 (Inada 
et al. 2008). Being able to distinguish between the two 
cases is important, not only for the understanding of the 
SDSS J2222+2745 lens system itself, but also in the con- 
text of using the statistics of cluster-lensed quasars as a 
cosmological tool. Naked cusps, i.e., lensing configura- 
tions where the source lies on a tangential caustic that ex- 
tends outside a radial caustic, are extremely rare among 
galaxy-scale lenses, the only known candidate among 
~ 100 such lens systems being APM 08279+5255 (Lewis 
et al. 2002). However, due to their shallower central mat- 



ter density slopes, naked-cusp lenses are expected to be 
much more common for clusters than for galaxies. Oguri 
& Keeton (2004) predicted that 30 - 60% of all cluster 
lenses will be of this type, the fraction depending on the 
central slope and triaxiality of clusters. 

To further investigate the lensing system configuration 
of SDSS J2222+2745 we searched our initial MOSCA 
imaging for potential stellar sources with blue colors con- 
sistent with the quasar. A blue source, labelled "D" 
in Figure [T] was identified just south of, and partially 
superposea on, galaxy G2. The source was apparently 
stellar, but the modest exposure time and overlap with 
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Table 2 

SDSS DR8 astrometry, photometry, and photometric redshift estimates for the three brightest quasar images and three brightest galaxies 

in SDSS J2222+2745 



Object 


RA 

(h m s) 


Dec 
(d m s) 


u 

SDSS 


g 

SDSS 


r 

SDSS 


i 

SDSS 


z 

SDSS 


Photometric 
redshift 





A 22:22:09.04 +27:45:37.9 24.55+0.93 21.10±0.03 20.67+0.04 20.87+0.07 21.00+0.27 

B 22:22:08.80 +27:45:40.0 22.07+0.19 21.43+0.04 21.03+0.05 21.27+0.10 21.08+0.29 

C 22:22:07.90 +27:45:37.8 24.50+0.93 21.51+0.05 21.14+0.05 21.28+0.10 22.14+0.63 

Gl 22:22:08.57 +27:45:35.6 25.42+2.41 20.92+0.08 19.00+0.03 18.05+0.02 17.65+0.05 0.485+0.015 

G2 22:22:08.70 +27:45:33.0 24.33+1.01 22.13+0.08 20.86+0.05 19.88+0.04 19.33+0.07 0.528+0.157 

G3 22:22:08.62 +27:45:29.9 24.39+2.32 22.28+0.20 20.43+0.07 19.48+0.05 19.13+0.16 0.509+0.031 



Note. — Epoch: 2009 October 18.19 UT 



Image A 
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Figure 2. NOT/ALFOSC spectra of the three brightest quasar 
images A-C. The most prominent quasar emission lines, redshifted 
to z = 2.82, are indicated in the figure. 

G2 made this conclusion somewhat uncertain. A set of 
deeper ugri images in better seeing was obtained at the 
NOT in September 2012 (see Table [lj. These images 
revealed image D to have a clearly stellar point-spread- 
function (PSF), motivating further spectroscopic follow- 
up of this image. On 2012 September 16-17 UT, a series 
of 8 x 2400s integrations were made using ALFOSC, using 
the same setup as for the earlier observations described 
in § 13. 1| The telescope pointing was dithered along the 
slit between exposures to improve fringe suppression at 
the red end of the spectrum. The spectra were flux cal- 
ibrated based on observations of the spectrophotometric 
standard star BD+33 2642. 

The slit orientation is illustrated in Figure [3j providing 
simultaneous spectroscopic coverage of images B and D, 
galaxies G2 and G3, and the giant arc. Inspection of the 
reduced 2D spectrum (see Figure |4| reveals the strong 
spectral line in the quasar B image around 4640 A, cor- 
responding to Lyman-alpha at the redshift of the quasar. 
A Lyman-alpha line is also clearly evident at this wave- 
length in the spectrum of the D image, confirming its 
identification as a gravitationally lensed image of the 
quasar. More surprisingly, a Lyman-alpha line at the 
same redshift is also seen superposed on the spectrum 
of galaxy G3, providing clear indication of a fifth (E) 
quasar image spatially coincident with this galaxy. The 
extracted ID spectra of image D and E in a ~ 1200 A wide 
region surrounding the Lyman-alpha line of the quasar 
are plotted in Figure [5] A clear peak is seen at the lo- 
cation of the Lyman-alpha line in both the D and E im- 
age spectra. Redward of this line, there is some visible 
discrepancy in the slopes of the two spectra. This dis- 
crepancy is caused by the spatial superposition of image 



E and the red galaxy G3, the latter having a continuum 
which is rising towards longer wavelengths. 

Figure[6]shows a scaled subtraction of the second epoch 
i-band image from the g-band image (see Table [lj. The 
PSF of the i-band image was matched to that of the 
g-band image by convolving it with an elliptical gaus- 
sian kernel. The z-band image was then scaled so that 
the subtraction optimally removes the light of galaxies 
G1-G3. Any remaining flux, shown in the negative in 
Figure [6j is bluer than galaxies G1-G3, and should be 
considered candidate images of the quasar. 

In Figure |6j image D is obvious. A comparison of 
this residual image to the residual quasar images A-C 
shows that it has a PSF consistent with a point source, 
as expected. Additionally there is excess blue flux at the 
spectroscopically indicated position of image E; this flux 
is obviously extended however. Figure [7J which shows 
the it-band image with contours from the g-i subtraction 
overlayed, shows the reason for this; the excess blue light 
is a combination of light from the core of galaxy G3, and 
offset from this, excess blue light attributable to image 
E. 

Finally, in Figure [6] there is also an additional blue 
excess source, near galaxy Gl, which we label F in Fig- 
ures[6]and [7] This source is consistent with being a point 
source. One of the slit placements from the initial spec- 
troscopic observations includes this possible quasar im- 
age, and Figure [5] includes a spectrum extracted at that 
location. The extracted spectrum suggests the presence 
of Lya emission; taken together, the imaging and spec- 
troscopic data suggest that this source is a sixth image 
of the quasar. 

3.3. Photometry of Quasar Images 

To achieve the most accurate photometry of images of 
the quasar, we have constructed a model of the lensing 
cluster core using GALFIT (Peng et al. 2010). The in- 
put PSF for each of the models discussed below was in all 
cases constructed by fitting bright, isolated and unsatu- 
rated stars using a multi-component Moffat profile using 
GALFIT. Typically 3 Moffat components were used; we 
have explored the effects of PSF reference star selection 
and find that it introduces an uncertainty of less than 
0.01 magnitudes for the brighter quasar images, with no 
significant biases. 

The cluster core model was constructed initially using a 
stack of the second epoch g, r and i images. Sources were 
placed initially at locations of obvious galaxies, and at 
the six quasar image positions discussed above. Quasar 
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Figure 3. Position of the slit for spectroscopic observations of Une. The wavelength of Lyman-alpha at the z = 2.82 redshift of 

SDSS J2222+2745 using ALFOSC in September 2012. The field of the quasar is indicated by the vertical dotted line. To dampen 

view and orientation is the same as for Figure [T] The background small-scale noise the spectra have been smoothed with a Gaussian 

image is a lh g-band MOSCA image combining all exposures ob- filter of FWHM= 18A 
tained in this passband in September 2012 (see Table III. 
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A [A] 

Figure 4. Reduced 2D spectrum with slit orientation as shown in 
Figure [3] The brightest spectrum near the bottom is the B image 
of the quasar. The locations of other sources along the slit are 
indicated along the left margin. The D and E images (the latter 
superposed on the spectrum of galaxy G3) display clear Lyman- 
alpha emission lines at the same redshift as the B image. The 
image has been gaussian smoothed on a scale of FWHM= 12A in 
the spectral direction and FWHM= 0.4" in the spatial direction. 

images are PSF components, and all galaxies are Sersic 
profile components. We continued adding components to 
the model iteratively at the locations of any significant 
residuals revealed by the modeling process, until a final 
complete model incorporating all significant sources was 
achieved. This model was then translated to individ- 
ual epoch and filter images by transforming the image 
coordinates to match the image to be analyzed, fixing 
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Figure 6. Residuals after subtracting the second-epoch stacked 
i-band image from the g-band image. Prior to the subtraction, the 
j-band image was PSF-matched by convolution with an elliptical 
gaussian kernel and then scaled to match the g-band flux of galax- 
ies G1-G3. The field of view and orientation is the same as for 
Figure [T] To indicate the locations of the galaxies we also show 
intensity contours from the i-band image. Objects with excess 
blue light in the cluster core are candidate quasar images. Images 
A-E are clearly confirmed by spectroscopy and readily visible in 
this image subtraction; a possible sixth quasar image (F) is also 
indicated. 

those coordinates, and then refitting the model compo- 
nent magnitudes, sizes and shapes to the image to be 
analyzed. 

The modeling process described above reveals that the 
images of the quasar include a contribution from the un- 
derlying host galaxy, as is seen in the two other quasars 
strongly lensed by galaxy clusters (Inada et al. 2005; 
Oguri et al. 2012). The data show that the host galaxy 
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Figure 7. A zoomed view - 4x finer than Figure[6]- of the u-band 
image. Red contours, as in Figure|6] show the second epoch i-band 
image. Blue contours show the residual g-i shown in greyscale in 
Figure [6] Note that at the position of image E, the residual blue 
light is extended, and that it partially coincides with an excess of 
■tt-band light that is most prominent at the center of galaxy G3. 
The peak of the excess blue light in g-i is offset from the center of 
G3; we interpreted this light as the sum of an image of the quasar 
(E) and a blue core in galaxy G3. 




Figure 8. Left Panel: The residuals in a 40x40 pixel region from 
the g-band image resulting from a single PSF GALFIT model cen- 
tered on image A. Right Panel: the same region fit by a PSF+Sersic 
profile. Note that the PSF-only fit clearly indicates the presence 
of an underlying host galaxy component, and that it is oriented 
tangentially to the cluster center, as expected, and that this com- 
ponent is well fit by a single Sersic profile. 

is blue - unlike the red hosts seen in the other two quasar- 
cluster lenses - and is most clearly seen in the brightest 
(and thus most likely the most magnified) image A. Fig- 
ure 8 shows residuals from a GALFIT model for image 
A in which only a PSF component is used, and a second 
model consisting of a PSF component plus a single Sersic 
component; this second model is a much better fit. 

In the g-band, including the host galaxy makes the 
fitted quasar 0.086 magnitudes fainter. In our judge- 
ment, none of the other images of the quasar allow for 
clear modeling of the host galaxy component - image B 
is bright enough, but complicated by a second galaxy 
less than two arcseconds to the north, and images C- 
F are insufficiently magnified. Thus, for images B-D, 
we proceed by measuring the quasar brightness using a 
simple PSF-only fit in GALFIT, and then correct each 
measured magnitude fainter by the host galaxy contri- 
bution as measured from image A. Though the quasar 
and host magnifications are not expected to be identical 



across the images, and the quasar/host ratio will not be 
completely fixed due to temporal variation of the quasar 
(see below), this simple procedure nevertheless avoids 
systematic noise that would be introduced by attempting 
to model the host galaxy independently in multiple im- 
ages. For reference, in addition to the g-band above, the 
correction for the host galaxy in r and i is measured to 
be 0.084 and 0.071 magnitudes, respectively. Within the 
uncertainties, this correction is the same across bands - 
i.e., the host galaxy has a color comparable to the quasar. 
For the it-band the image is of insufficient depth to model 
the host galaxy directly, we adopt the same correction as 
the g-band. The resulting photometry for images A-D, 
across four bands and two epochs, is reported in Table 3. 
Uncertainties are calculated by using the full GALFIT 
model to generate realisations of each measured frame 
with appropriate noise, which are then refit. The mea- 
sured scatter in the photometry of quasar images, com- 
bined in quadrature with the estimated calibration un- 
certainty is reported. For the brighter images, system- 
atics due to uncertain treatment of the host galaxy are 
likely dominant, and not addressable with these limited 
data. 

Images E and F must be treated with further care. The 
photometry of image E is complicated by the blue core 
of galaxy G3, and image F is simply very faint, and near 
the core of the brightest cluster galaxy, Gl. Particu- 
larly at redder wavelengths, this makes its measurement 
highly uncertain. Furthermore, tests using the full GAL- 
FIT model, in which we re-measured noisy realisations 
of the full model, also show potentially strong biases - 
up to factors of 3 - in the photometry of image F. To at- 
tempt best possible photometry for these we proceeded 
as follows. 

First, note that if the colors of the various quasar 
images are the same, then the relative photometry of 
residual images in the g-i subtraction in Figure [6] can 
be used to measure the fainter images. Matched colors 
are not guaranteed, due both to intrinsic variability (e.g., 
Schmidt et al. 2012) coupled with time delays, as well as 
variations induced by differential magnification between 
the quasar and underlying host galaxy, and possible ex- 
tinction of individual images due to intervening galaxies 
in the complex cluster core. However, we can validate the 
likely variability in color between the images by consider- 
ing the three brightest and most isolated images (A-C); 
these show that the typical scatter about a fixed color 
is 0.05 in g-i and 0.03 in g-r. This uncertainty is much 
less than the systemic biases noted above in attempts at 
direct measurement, and so the fixed color assumption is 
sufficient given these data. 

Thus for image F, we use difference images to mea- 
sure the differential flux of image F relative to images 
A-C. Assuming that F has the same color as these im- 
ages, the flux of image F is then given simply by the 
mean flux of images A-C times the flux ratio in the dif- 
ference image. All magnitudes are measured using PSF 
fitting photometry, for consistency. The difference image 
g-i provides the best measurement of the residual fluxes, 
and is the basis of the measurements reported in Table 3. 
As before, we estimate uncertainties by using the fitted 
GALFIT model to generate and refit noisy realisations 
of the measured image. In the case of image F, the final 
reported uncertainty is the scatter in the magnitude off- 
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set in the residual image, combined in quadrature with 
the color scatter between the brighter images. 

For image E we must further account for the residual 
flux from the core of galaxy G3. We proceed as for image 
F, but additionally use GALFIT to model a best-fitting 
Sersic profile fixed at the center of G3 concurrent with 
the PSF component for image E. 

Finally, we note that each frame is directly calibrated 
to the SDSS using many hundreds of objects; calibration 
uncertainties are 0.01 magnitudes at worst. 

3.4. Lensing cluster 

The maximum separation of the quasar images is 15'.'1 
(between quasar images A and C) clearly indicative of a 
cluster-sized lensing mass. The SDSS photometric red- 
shift estimates of galaxies G1-G3 (see Table [2]) are con- 
sistent with all three galaxies being located m a clu ster 
at z ~ 0.5 . As part of the observations described in § |3.1| 
and § |3.2| above, all three galaxies were targeted spec- 
troscopically. Redshift estimates were obtained through 
cross-correlation with the elliptical galaxy template of 
Kinney et al. (1996), yielding the redshift estimates listed 
in Table SJ The redshifts are consistent within the mea- 
sured uncertainties, and all three galaxies G1-G3 clearly 
belong to the same cluster at a redshift consistent with 
the SDSS photometric redshift estimates. To further in- 
vestigate the properties of the lensing cluster, we applied 
the SExtractor software (Bertin & Arnouts 1996) for 
galaxy photometry to the stacked second-epoch MOSCA 
gri images. Prior to obtaining the photometric measure- 
ments, the stacked r- and i-band images were smoothed 
to match the resolution of the <?-band image. Object de- 
tection was initially performed on the stacked g-band 
image, and SExtractor was subsequently run in dual- 
image mode on the stacked images in the other pass- 
bands, using the g-band image as a reference for object 
detection. The observed magnitudes were subsequently 
converted to AB magnitudes in the SDSS system using 
objects in the MOSCA field with catalogued SDSS pho- 
tometry. Total magnitudes were determined using the 
SExtractor MAG_AUT0 parameter, and colors were deter- 
mined from aperture magnitudes (MAG_APER) within a 2" 
diameter aperture. A g-i vs. i color-magnitude diagram 
of non-stellar objects, selected based on their location in 
a magnitude - /x max diagram (where /x max is the central 
surface brightness of objects), is shown in Figure [9] A 
red sequence is evident at g-i ~ 2.5, similar to the g-i 
colors of the three galaxies G1-G3 (see Table 4]). Fig- 
ure [TO] shows a smoothed map of the light distribution of 
red sequence galaxies over the observed MOSCA field. A 
strong peak in the luminosity distribution is located at 
the position of galaxies G1-G3. 

3.5. Giant arc 

The identification of SDSS J2222+2745 as a candidate 
gravitational lens system was partially due to recognition 
of the giant arc just south of the cluster center. The 
extent of the arc in our deepest stacked g-band image (see 



Figure |3j) is 10". Our SExtractor photometry (see § 3.4 1 
of the arc yields g = 21.57, and u-g — 0.87, g-r — 0.21, 
g-i = 0.41, where the colors are measured in 2" diameter 
apertures centered on the brightest central portion of the 
arc. We note that the quoted magnitude of the arc can be 
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Figure 9. Color-magnitude diagram of 3033 non-stellar objects 
detected in the MOSCA field containing SDSS J2222+2745. The 
galaxies Gl-3 are indicated by green dots. Other galaxies with g-i 
colors within ±0.35 magnitudes from the red sequence formed by 
early-type galaxies in the lensing cluster are indicated by red dots. 
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Figure 10. A contour plot showing a smoothed map of the light 
distribution of the red sequence galaxies indicated by the red and 
green colors in Figure [9] The background figure is a gri composite 
based on the second-epoch MOSCA images. The coordinates are 
relative to the position of galaxy Gl, which closely coincides with 
the cluster lens (see § |4j . A strong peak in the light distribution is 
seen at this location. 

considered a lower limit on its brightness, given the very 
extended mophology of the arc. Although SEXtractor 
attempts to treat the arc and the WD as separate objects, 
a more careful modeling of the shape of the arc (including 
PSF modeling and removal of the superposed white dwarf 
point source) would be required to obtain an accurate 
total magnitude. The spectrum of the arc shows a single 
emission line at A b s =~ 4015A, which we interpret as 
Lyman-alpha at z = 2.30. The alternative interpretation 
of [O II] (3727 A) at z = 0.077 would put the arc in front 
of the cluster, which is obviously inconsistent with the 
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Table 3 

Photometry of QSO images 



Data Frame Quasar Image 

A B C D 



(/-band 1st epoch 
cj-band 2nd epoch 


21.03±0.01 
21.15±0.01 


21.52±0.01 
21.57±0.01 


22.13±0.01 
22.22±0.01 


23.48±0.04 
23.74±0.03 


24.32±0.09 
24.22±0.08 


25.21±0.24 
24.81±0.07 


r-band 1st epoch 
r-band 2nd epoch 


20.66±0.01 
20.80±0.01 


21.13±0.01 
21.16±0.01 


21.69±0.01 
21.79±0.01 


23.13±0.11 
23.46±0.08 


23.92±0.09 
23.82±0.08 


24.81±0.24 
24.41±0.07 


j-band 1st epoch 
j-band 2nd epoch 


20.82±0.01 
20.95±0.01 


21.20±0.01 
21.31±0.01 


21.81±0.01 
21.90±0.01 


23.46±0.12 
23.78±0.08 


24.04±0.09 
23.96± 0.08 


24.92±0.24 
24.55±0.07 


M-band 2nd epoch 


22.85±0.08 


23.29±0.09 


24.18±0.17 


>24.58 


>24.58 


>24.58 



Uncertainties are quoted at the 1-cr level. Images D,E,F are not detected in the M-band; the 1-a detection limit in our data is 24.58 mag. 
When comparing this photometry to the SDSS photometry reported in Table 2, a correction for the quasar host galaxy should be applied 
to the SDSS photometry, as detailed in the main text. 



Table 4 

Photometry and redshift measurements for galaxies in the cluster core. 



Galaxy 


i 


g-i 

(2'.'0 aperture) 


r— i 

(2'.'0 aperture) 


Redshift 




Gl 
G2 
G3 


18.74 
18.97 
19.37 


2.71 
2.53 
2.45 


0.97 
0.96 
0.93 


0.493 ±0.013 
0.496 ±0.012 
0.489 ±0.013 



lensed morphology of the arc. The lens model presented 
in the following section does not predict a counter-image 
of the arc. 

4. LENS MODEL 

We compute a parametric lens model using the publicly 
available software Lenstool (Julio et al. 2007), through 
Markov Chain Monte Carlo minimization in the image 
plane. The lens is modeled with several components: a 
cluster halo, BCG, and cluster galaxies. All the com- 
ponents are represented by a pseudo-isothermal ellipti- 
cal mass distributions (PIEMEp"[ Limosin et al. 2005, 
Ehasdottir et al. 2007), parametrized by its position x, 
y; a fiducial velocity dispersion a; a core radius r core ; a 
scale radius r s ; ellipticity e = (a 2 — b 2 )/(a 2 + b 2 ), where 
a and b are the semi major and semi minor axes, respec- 
tively; and a position angle 9. All the parameters of the 
cluster halo distribution are allowed to vary within pri- 
ors, except for r s , which is not constrained by the lensed 
images. 

Cluster members are selected based on their g-i color 
in a color-magnitude diagram, with respect to the cluster 
red-sequence (see FigureM and their position, ellipticity 
and position angle are fixed at their observed values, as 
measured with SExtractor. The parameters a and r s are 
determined through scaling relations (see Limousin et 
al. 2005 for a description of the scaling relations). Three 
galaxies, Gl, G2 and G3, are in close proximity to three 
lensed images of the quasar, and thus affect their po- 
sitions more directly than other cluster members. To 
account for the intrinsic scatter in the scaling relations, 
we allow some of the parameters of these three galax- 
ies to be modeled separately from the rest of the cluster 
galaxies. 

11 This profile is formally the same as dual Pseudo Isothermal 
Elliptical Mass Distribution (dPIE, see Eh'asdottir et al. 2007). 



As constraints, we use the redshift and positions of the 
images of the quasar and of the lensed galaxy at z = 2.30. 
We have a total of 12 free parameters: 6 free parameters 
for the cluster halo, and 2 free parameters for each of 
the galaxies, and 14 constraints. The best fit parameters 
are listed in Table [5] The uncertainties were estimated 
through the MCMC sampling of the parameter space. 
For parameters that were not well-constrained, we list 
the priors that were set in the minimization. To account 
for systematics due to misidentification of the sixth im- 
age, we also compute a model assuming only five images. 
The two families of models are similar, although we note 
that Gl is poorly constrained in the model without image 
F. 

The image-plane RMS scatter of the quasar images is 
0.14" in the 5-image model and 0.18" in the 6-image 
model; the RMS of the blue arc is somewhat larger 
(0.26", 0.38", respectively), due to the difficulty of mea- 
suring its positional constraints, and the larger weight 
given to the quasar images in the modeling process. 

Both models successfully predict the locations of all 
the observed images (either five or six, depending on the 
model constraints), and do not predict additional im- 
ages, with the exception of demagnified central images 
that cannot be detected in our data: the model with 
six images predicts a demagnified seventh image close 
to the center of galaxy G2. This is expected in this 
strong lensing configuration. Our imaging data do not 
reveal any hint of this seventh image, though a positive 
identification might be possible with the superior spatial 
resolution of HST. A careful re-examination of the spec- 
troscopic data for G2 shown in Figure [4] and plotted in 
Figure [5] may indicate a weak signature of the Lyman- 
alpha line from the seventh image (G), but this feature is 
almost at the level of the noise, and more sensitive spec- 
troscopy with 8-10m class telescopes will be required for 
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Figure 11. Caustics, critical lines, and image positions from 
the best-fit model, superposed on a gri color image of the 
SDSS J2222+2745 lens system. The critical curves for a source 
at the quasar redshift (z = 2.82) are indicated in red. The caustic 
curves in the source plane are shown as yellow lines, and the pre- 
dicted source locations of the quasar and arc with respect to these 
curves are indicated by green crosses. The observed locations of 
the six quasar images A-F are indicated by blue circles. 

a definite detection. 

Table [6] lists the model-predicted magnifications and 
time delays for the two model families. Uncertainties 
are computed from a suite of models, drawn from the 
MCMC chain, that represent the 1-ct uncertainty in the 
parameter space. All time delays are given relative to 
image A. The best-fit models predict short time delays 
of order of 100 days between the close pair (A-B), and 
between the central images (D-E-F) . Image C precedes A 
and B by 1400 days, and D-E-F follow A and B by 700 - 
1000 days. In the 6-image model, 66% of the models, as 
well as the best-fit model, predict the order of the images 
to be C-B-A-F-E-D. All the models predict that image C 
is the first one to occur, followed by B-A (in 93% of the 
models) or A-B (7%); followed by F-E-D (72%) F-D-E 
(16%) or D-F-E (12%). In the 5-image model, 62% of 
the models, as well as the best-fit model, predict C-B-A- 
E-D; here too, all the models predict that C is the first 
image to occur, followed by B-A (77%) or A-B (23%), 
and then E-D (81%) or D-E (19%). 

5. SUMMARY AND DISCUSSION 

5.1. Previous observations and predictions of high 
image multiplicities 

We have identified SDSS J2222+2745 as a gravita- 
tional lens system producing six observed images of the 
same quasar. This is the third known cluster-scale lens 
and the only such lens producing more than five observed 
images. The only previously known six-image AGN lens 
is B1359+154 (Rusin et al. 2001), where a radio source 
at z s = 3.235 is being lensed by a compact group of 
three galaxies at zi ~ 1. However, the quasar image sep- 



arations we observe in SDSS J2222+2745 are an order of 
magnitude larger than in B1359+154, and the lensed im- 
ages are ~ 3 magnitudes brighter at optical wavelengths, 
making six images detectable in ground-based imaging of 
SDSS J2222+2745. The B1359+154 lensed image config- 
uration is similar to what we find in SDSS J2222+2745, 
with the three fainter images located in close proximity 
to the three galaxies situated within the Einstein radius 
of the lens system, and three brighter images located fur- 
ther away from these galaxies. However, the lens model 
of B1359+154 presented by Rusin et al. (2001) does not 
allow a significant dark matter halo associated with the 
galaxy group in addition to the three halos corresponding 
to the three individual galaxies. In contrast, in our model 
the cluster halo component is the dominant contributor 
to the observed image splitting of SDSS J2222+2745. 

Since cluster-scale quasar lenses and image multiplici- 
ties > 5 are both rare, we would expect a lensing system 
such as SDSS J2222+2745 to be extremely rare. There 
are no appropriate quantitative predictions for the oc- 
currence of this type of lens systems available in the lit- 
erature. Evans & Witt (2001) calculate the incidence 
of sextuplet and octuplet image systems (with demagni- 
fied seventh and ninth images, respectively) and estimate 
that these constitute ~ 1% of the population of quasars 
multiply imaged by galaxy-scale lenses. This is consis- 
tent with the single known sextuplet case of B1359+154 
among the ~ 100 known galaxy-scale quasar lenses. For 
the larger image separations produced by cluster lenses, 
published studies of the probabilities of different lensing 
multiplicities (Oguri & Keeton 2004; Li et al. 2007; Mi- 
nor & Kaplinghat 2008) only include double, quadruple 
(with demagnified third and fifth images, respectively) 
and naked cusp configurations. Li et al. (2007) state 
that they do not expect a significant fraction of sources 
to have image numbers exceeding 5. 

5.2. Comparing the lens and source redshifts to 
predictions 

Both the quasar source redshift and cluster lens 
redshift of SDSS J2222+2745 are comparable to the 
two other cluster quasar lenses and in good agree- 
ment with theoretical expectations: The observed lens 
redshift in SDSS J2222+2745 is z = 0.49, which is 
similar, but slightly lower than the cluster lenses in 
SDSSJ1004+4112 \z = 0.68) and SDSS J1029+2623 
[z = 0.58). This redshift is very close to the peak (at 
z ~ 0.5) of the predicted lens redshift distribution for 
photometrically selected wide-separation quasars in the 
SDSS (Li et al. 2007). The observed quasar source red- 
shift in SDSSJ2222+2745 is z = 2.82, slightly higher 
than the sources in SDSS J1004+4112 {z = 1.73) and 
SDSS J1029+2623 (z = 2.20). These three redshifts are 
all near the peak (around z = 2) of the predicted source 
redshift distribution for photometrically selected wide- 
separation quasars in the SDSS (Li et al. 2007). 

5.3. Prospects for time delay measurements 

The photometry reported in § |3.3| shows significant 
variability in the quasar, demonstrating the feasibility of 
making time delay measurements in SDSS J2222+2745 
from monitoring campaigns using ground-based 2-4m 
class telescopes. We have recently initiated a monitoring 
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Table 5 

Best-fit lens model parameters 



Halo 
(PIEMD) 


RA 

(") 


Dec 
(") 


e 


e 

(deg) 


^corc 

(kpc) 


»"cut 

(kpc) 


cro 

(km s ) 








Model 1: 


constraints from 5 QSO images. 


RMS QSO =0.14" 






Cluster 
Gl 
G2 
G3 

L* galaxy 


n i * +0.50 
-°- 14 -1.03 

[0.0] 
[-0.596] 
[-1.762] 


-, ,o+1.90 
1-1S-0.81 

[0.0] 
[-5.220] 
[-2.308] 


fi6+ 0,25 
[0.26] 
[0.21] 
[0.41] 


97+5 
[3.8] 
[75.5] 
[-18.3] 


Ofi+15 
26_ 42 

77+ - 77 
u -' '-0.23 

0.13 [0.0 0.3] 
0.19 [0.0 0.3] 
[0.15] 


[ZUUUJ 

[47] 
[24] 
[20] 
[20] 


eon+64 
&8U -163 

155+ 55 
100 _168 
91 9+40 

279+5? 

[100] 






Model 2: constraints from 6 QSO images. 


RMS QSO =0.18" 






Cluster 
Gl 
G2 
G3 

L* galaxy 


-0 16+ 058 

U ' ±D -0.71 

[0.0] 
[-0.596] 
[-1.762] 


2-10l -.90 
[0.0] 

[-5.220] 

[-2.308] 


75+ 013 
U - ' S -0.05 

[0.26] 
[0.21] 
[0.41] 


96 + 3 
[3.8] 
[75.5] 
[-18.3] 


Q7+13 
°' -40 
05+ - 04 

0.15 [0.0 0.3] 

0.15 [0.0 0.3] 

[0.15] 


[2000] 

[47] 
[21] 
[20] 
[20] 


578±f 14 
303^ 
226+| 

2 35_ 33 

[100] 



Note. — All coordinates are measured in arcseconds relative to the center of Gl, at [RA, Dec] = [335. 53573, 27.75986]. The ellipticity 
is expressed as e = (a 2 — b 2 )/(a 2 + b 2 ). 9 is measured north of West. Error bars correspond to 1 a confidence level as inferred from the 
MCMC optimization. Values in square brackets are for parameters that are not optimized. Errors in square brackets represent the lower 
and upper limits that were set as prior in the optimization process, for parameters that were not well-constrained by the data. The location 
and the ellipticity of the matter clumps associated with the cluster galaxies were kept fixed according to their light distribution. 



Table 6 







Model-predicted ma; 


mifications and time delays 




Image 


Magnification (1) 


Time delay (1) [days] 


Magnification (2) 


Time delay (2) [days] 


A 
B 
C 
D 
E 
F 


c 7 +14.5 
°-'-3.1 

5-6t 2 2 2 8 

9 c+3.6 

6+ - 7 
U.D_ 2 

7+0.9 




s7 +187 
°'-296 
-1399177° 

960±1<!5 

ocri +448 
° oi -278 


4.71?;* 

5-4lt' 6 6 
9 4 +l.l 
Z - 4 -0.3 
7+O.3 
u -'-0.2 

6+°' 3 

U.D_q 2 

5+°' 5 
U.O_ 2 




1,9 + 158 
±lz -225 

-1368l«6 

9311233 
806li° 4 
723l 2 g 3 



Note. — Model (1): constraints from 5 QSO images. Model (2): constraints from 1 
to image A. 



QSO images. Time delay is given in days, relative 



program for SDSS J2222+2745 at the NOT for this pur- 
pose. The predicted time delays from our lens model (see 
§ [I| indicate that three independent time delays AtAB 
(= AtB — At a), At ed and AtFE in the system are all 
of order 100 days. This suggests that several time de- 
lay measurements could, in the best case, be obtained 
from a single season of monitoring. The predicted time 
delays between image C and the close pair (A-B) and 
between (A-B) and the central images (D-E-F) are of 
order 1400 days and 700-1000 days, respectively. These 
values are similar to the two longest predicted time de- 
lays for SDSS J1004+4112, of 1218 and 1674 days, re- 
spectively (Oguri 2010) and the measured longest time 
delay in SDSS J1029+2623 of At AB = (744 ± 10) days 
(Fohlmeister et al. 2012). 

5.4. Future work 

Morphological, photometric and spectroscopic evi- 
dence presented in this paper strongly suggest the pres- 
ence of a sixth (F) image of the quasar located closest 
to the lensing mass center. The odd number theorem of 
strong gravitational lensing theory will then require the 



presence of a seventh (G) image, for which our model 
predicts a location close to the center of cluster galaxy 
G2. In the case of galaxy-scale lenses, the faintest image 
of a multiply lensed quasar is typically highly demagni- 
fied and hence unobservable in the vast majority of cases. 
However, the shallower central slopes of the mass density 
profiles of clusters increases the chance of observing the 
final odd image, as exemplified by the identification of 
the demagnified fifth central image in SDSS J1004+4112 
using HST imaging (Inada et al. 2005, 2008). 

While the data presented in this paper do not provide 
conclusive evidence for the presence of the seventh im- 
age, there is a hint of Lyman-alpha emission from this 
image in the spectrum of G2. If this feature is indeed 
emission from the quasar, a firm confirmation would be 
possible through spectroscopic observations with 8-10m 
class telescopes. 

While time delay measurements would provide 
important additional modelling constraints on the 
SDSS J2222+2745 lens system, follow-up HST imaging 
will also be crucial to improve the modeling constraints 
from the observed locations of the quasar images and 
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hopefully provide morphological and photometric confir- 
mation of the predicted seventh (G) image close to the 
center of galaxy G2. Further constraints on the lens sys- 
tem could also come from spatially resolved, multiply im- 
aged structures within the giant arc and from the identifi- 
cation of additional fainter, multiply imaged, galaxies as 
revealed by the HST observations of SDSS J1004+4112 
(Sharon et al. 2005) and SDSS J1029+2623 (Oguri et al. 
2012). 

5.5. Broader significance of this discovery 

The significance of the discovery of SDSS J2222+2745 
is twofold: Firstly, the presence of six quasar images 
(and potential future confirmation of a seventh image), 
spread over a range of clustercentric radii, with good 
prospects for future time delay measurements, coupled 
with the presence of an additional lensed source pro- 
ducing a bright giant arc, can provide unique detailed 
measurements of the matter distribution in the core of a 
cluster. Such measurements provide a stringent test of 
the scenario of structure formation, including the prop- 
erties of dark matter in cluster cores. 

Secondly, adding our new discovery of 
SDSS J2222+2745 to the previous detections of 
SDSSJ1004+4112 and SDSS J1029+2623 makes a 
significant contribution towards establishing quasars 
strongly lensed by galaxy clusters as a cosmological 
tool. The probability of such rare lenses depends very 
sensitively on cosmological parameters, in particular 
the normalization of the matter power spectrum 
and the matter density f2 M (e.g., Li et al. 2007). Even 
though no directly applicable predictions can be found 
in the literature, the occurrence of lens systems with 
image separations and image multiplicities as high or 
higher than SDSS J2222+2745 must necessarily be 
exceptionally rare, and this could conceivably be the 
the only such system, provided that the current ACDM 
paradigm is correct. 
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